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The two-dimensional electron gas at the LaTiOa/SrTiOs or LaAlOa/SrTiO;} oxide interfaces be- 
comes superconducting when the carrier density is tuned by gating. The measured resistance and 
superfluid density reveal an inhomogeneous superconductivity resulting from percolation of fila- 
mentary structures of superconducting "puddles" with randomly distributed critical temperatures, 
embedded in a non-superconducting matrix. Following the evidence that superconductivity is re- 
lated to the appearance of high-mobility carriers, we model intra-puddle superconductivity by a 
multi-band system within a weak coupling BCS scheme. The microscopic parameters, extracted by 
fitting the transport data with a percolative model, yield a consistent description of the dependence 
of the average intra-puddle critical temperature and superfluid density on the carrier density. 



PACS numbers: 74.81.-g,73.20.-r,73.40.-c,74.20.- 



The discovery of superconductivity in the two-dimensional electron gas (2DEG) formed at certain oxide interfaces 
such as LaA10 3 /SrTi0 3 (LAO/STO) or LaTi0 3 /SrTi0 3 (LTO/STO) pJS] has attracted great interest and stimulated 
intense research activity. The possibility of tuning the carrier density by means of electrostatic gating in these 
extremely two-dimensional superconductors (the Fermi wavelength being smaller than the thickness of the gas [3]) 
opens new opportunities to study fundamental issues in quantum fluids [5]. More generally, these 2DEG appear as 
very rich quantum matter systems, since they also present magnetism [BH5] and strong tunable spin-orbit coupling 
[113 E]- As for every two dimensional systems, disorder is expected to play a relevant role on their physical properties. 
Recent magnetotransport measurements revealed the existence of two kinds of carriers in LTO/STO, with high and 
low mobility respectively [J- It has been suggested that superconductivity is related to the appearance of the high- 
mobility carriers[H 112] with an inhomogeneous character, i.e. superconducting "puddles" embedded in a (weakly 
localizing) metallic background [5 . Direct superfluid density [5] and magnetic [5] measurements display signs of strong 
inhomogeneity at micrometric scales. Recent current distribution^ and surface potential studiesfjjl] reveal stripes 
structures in these interfaces. Given all these signs of inhomogeneities, one has first to analyze the role of disorder 
before getting further insight in the physics of the 2DEG, and for example, explain the rapid decrease of T c with the 
gate voltage while under-doping. 

Our aim is on the one hand to infer a more detailed structure of the inhomogeneous state and on the other hand 
to provide a description for superconductivity in an inhomogeneous system with different kinds of carriers. We first 
show that resistance measurements [4j and the topographic mapping of the superfluid density on micrometric scale [9] 
can be phenomenologically accounted for within a scheme of a percolative inhomogeneous 2DEG with poor long- 
distance connectivity, but substantial superconducting (SC) current loops to warrant a high fraction of diamagnetic 
response. We then show that the properties of the SC puddles (e.g., their fraction, and critical temperatures) can be 
extracted from experiments and used at a microscopic level to model the intra-puddle multi-carrier pairing, gaining 
insight about the pairing mechanism. Although some features of the diamagnetic response are seemingly related to 
a strong coupling SC behavior[9], our analysis demonstrates that inhomogeneities and multi-band superconductivity 
fully account for the SC behavior of these systems within a standard weak coupling BCS scheme. 

— Sheet resistance — We have recently shown that the percolative SC transition in inhomogeneous systems, like the 
oxide interfaces, is well described by the Effective Medium Theory (EMT) [15] • For samples exhibiting a resistance 
drop due to superconductivity, the measured sheet resistance of LTO/STO is well fitted by the EMT with a SC 
fraction w < I (the SC puddles). Each puddle has a random local critical temperature T c , which we assume to have 
a Gaussian distribution, with mean T c and width 7. The weight w of the T c distribution represents the fraction of 
SC puddles and the non SC fraction 1 — w represents the metallic background (we refer the reader to Ref. [T5J for 
details). The sheet resistance at any temperature is 



R(T) = R c 



'T — T c 

> erf ( — — — - ) + 1 — w 

V2j 



with fitting parameters R 00 (high-temperature sheet resistance), to, T c , and 7. The percolative SC transition, marked 
by a vanishing R, takes place at a temperature T p < T c , such that 

This equation has a solution for T p only if w > 5, i- e -i if the SC fraction can percolate in the two-dimensional 
system. When T p becomes negative or is not defined, the system remains resistive down to T = 0, although a 
significant reduction of sheet resistance may still occur, witnessing the presence of a sizable (though not percolating) 
SC fraction. 

Fig. [l] (a) reports the measured sheet resistance as a function of temperature for various values of the gate voltage 
V g in a LTO/STO sample, and the fitting EMT curves. From the fits one can extract the average intra-puddle critical 
temperature T c [corresponding to the temperature of maximal slope of R(T) within EMT], which is reported in Fig. 
[31(a) as a function of V g (red solid line and empty circles). The fits also provide the fraction of SC regions w [displayed 
in Fig. [I] (c) as a function of V g ] . Since EMT disregards spatial correlations, the simultaneous request of a percolating 
SC cluster and the tailish shape of the resistance near percolation force the weight of the regions which can become 
SC to be w « \ [15] (Fig. [l[c), see also Supplemental Material). Fig. |l[d) reports the width 7 of the Gaussian 
distribution of T c , which diverges as the fraction w of the SC puddles goes to zero. This is rather natural because the 
density decrease emphasizes the effects of disorder so that fluctuations of the local superconductivity increase, leading 
to a substantial broadening of the T c distribution. We will show later a more precise description of this phenomenon 
within the two-bands model. 
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FIG. 1. (a) Sheet resistance as a function of temperature in a LTO/STO sample for various values of the gate voltage. From 
top to bottom V g — —10, 10, 30, 50, 70, 100, 150, 200 V. Symbols mark the experimental data, while the dashed lines are the 
theoretical fits using the EMT (see text), (b) From top to bottom V g — 50, 70, 100, 150, 200 V. Symbols mark the experimental 
data, while the dashed lines are the theoretical fits using the RRN (see text). Inset: sketch of an inhomogeneous RRN system 
with a SC filamentary+baubles cluster. The shaded region is a closed loop contributing to the diamagnetic response, (c) SC 
fraction w extracted from the fits of the sheet resistance within EMT (filled circles, solid line) and within the RRN (empty 
diamonds), (d) Width of the Gaussian distribution of T c used in the fits of the sheet resistance within the EMT (filled circles, 
solid line) and within the RRN (empty diamonds). 



- Superfluid density — The superfiuid density was measured in LAO/STO interfaces using a scanning SQUID 
probe averaging over micrometric regions [5]. This technique is not sensitive to submicrometric inhomogeneities, but 
revealed a distribution of superfiuid densities within a given sample on the micrometric scale, thereby supporting the 
idea of inhomogeneous 2DEG at these oxide interfaces. Since the sheet resistance curves of these systems are similar 
to those of LTO/STO interfaces, we assume that our percolative EMT analysis also applies in this case and that the 
measured local superfiuid density is an average over an inhomogeneous state of submicrometric puddles |17j . 

To capture this effect, we extend the EMT to small but finite frequency w. We model the metallic regions with 
a Drude-like complex conductivity o~n(uj) = B/(A + iu>) and the SC regions with a purely reactive conductivity 
o~s(uj) — B/(iu>). We introduce the quantities ps(w) = \/as{uj) — ico/B and Pn{u) = 1/o~n(w) = po + ps(u>), with 
po = A/B. At high temperature, the system is metallic and p{oS) = 1/o~(uj) — pn everywhere. However, within the 
SC puddles, pg vanishes as soon as the local critical temperature is reached. The solution of the EMT equation, 
neglecting higher frequency terms ~ w 2 , is 



p(u) « Po(w N - ws)ti(w N - w s ) + 



Ps(u) 

\w N -w s \' 



where -d(x) is the Heaviside step function, ws is the fraction of puddles that have become SC (at a given temperature 
T), and iojv = 1 — u>s is the metallic fraction (resulting both from puddles that have not yet become SC and from the 
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FIG. 2. (a) EMT fits of the superfiuid density in LAO/STO (lines) as a function of temperature for various gate voltages. 
The experimental data (symbols) are extracted from Ref. 9, (b) SC fraction w responsible for the diamagnetic response, (c) ) 
Width of the Gaussian distribution of T c used in the fits of the superfiuid density. 



metallic background). Evidently, when wn > ws, the system displays a Drude-like complex conductivity. However, 
below the percolation temperature T p (if any), ws > wn, and we find a purely reactive conductivity 



a(u) 



A(w s - w N ) 



P(w) 



Therefore, for T < T p and w given by Eq. (TT|) , the superfiuid density of the percolating network is 



J s oc w s - w N 



1 — w erf 



erf 



T P -T C 



-erf 



T-T c 
T-T, 



(2) 
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Fig. [2] (a) reports this micrometrically averaged superfiuid density together with our EMT fits [Eq. p])] at various 
gate voltages V g . Remarkably, the T dependence of J S (T)/J S (T = 0) mimics the qualitative behavior of the BCS 
prediction, but may quantitatively deviate from it. The slope at the transition, for instance, is ruled by the width 7 
of the distribution of critical temperatures. We can thus account for the deviations from standard BCS observed in 
Ref. [3 which were attributed to a tendency to strong coupling. Our alternative explanation for these deviations is 
that the system is in a weak coupling regime, but the percolative superfiuid density is ruled by the submicrometric 
distribution of inhomogeneities averaged at the micron scale by the SQUID pick-up loop used in Ref SI The fits allow 
to obtain the fraction w of SC puddles extracted from diamagnetic measurements which is reported in Fig. ^lb). It 
ranges from 1/2 to 1, and is therefore always larger than the fraction obtained from transport measurements. This 
is fully consistent within a percolative description of the system: resistance measurements mainly probe the direct 
percolative path, regardless of dead-ends and disconnected SC regions[H], while screening measurements are sensitive 
to all SC loops, even when not connected to the backbone. In this case, the diamagnetic fraction can be nearly one, 
while the connectivity for macroscopic transport can be very small (or even vanishing) if many SC puddles or loops 
are disconnected [see the sketch in the inset of Fig. [lib)]. Fig. W(c) displays the behavior of the width 7 of the 
T c distribution. Despite the differences in material and physical quantities, 7 has a similar behavior to the 7 from 
transport in LTO/STO [Fig. [lid)] and is of comparable magnitude. 

Effect of space correlations — The explanation of the above discrepancy between the SC weight w in transport 
and magnetic experiments rests on the filamentary structure of the SC fraction in LXO/STO systems, which is one 
of the two main outcomes of our work. However, the above results (and fits) have been obtained within the EMT, 
which is a mean-field-likc approach neglecting space correlations. Therefore, to substantiate the above interpretation, 
it is important to test how robust is the finding that the SC fraction involved in transport is rather small, w ~ 0.5, 
within a model, which takes into account the space correlations inherent in the filamentary structure. To this 
purpose, we refine our model and solve a random resistor network (RRN) where the SC regions only occur on a 
cluster with strong space correlations embedded in a metallic matrix. Based on a preliminary analysis [IB], showing 
that the tailish sheet resistance implies both a dense SC cluster at short distances and a filamentary structure at 
larger distance, we numerically generate a fractal-like structure |181 ITU] . Due to the poor long-distance connectivity, 



£ k,C — £ £ + ~ — — + 



percolation only occurs when the resistance of almost all bonds has been switched off. For a Gaussian distribution of 
critical temperatures, the system is thus forced to explore the low temperature asymptotic of the distribution, thereby 
producing a tailish resistance [20]. On the other hand, a purely filamentary structure, although dense at short distance 
due to the almost fractal cluster, is still too faint (w < 0.3) and fails to reproduce accurately the high temperature 
part of the R(T) curves. Thus in order to vary the density of the SC cluster without crucially changing the long-range 
connectivity, we decorate the filaments with randomly distributed circular "superpuddles" , their number being chosen 
to yield RRNs with weights w ranging from 0.3 to 0.7 (see Supplemental Material). The resulting clusters look like 
a filamentary "Christmas tree" with decorating baubles [see inset of Fig. [TJ (b)]. Fig. [TJ (b) reports fits of the sheet 
resistance using resistivity curves calculated on such Christmas trees. Noticeably, even though the SC weight w is not 
forced to | like in EMT, we still find that the sheet resistance only acquires a tailish behavior if 0.5 < w < 0.65 (Fig. 
IT] (c)), where the lower bound is imposed by the high slope at intermediate temperatures and the upper bound is due 
to the tailish behavior close to percolation. Fig. [TJ (d) reports the width 7 of the Gaussian distribution of T c obtained 
within the Christmas tree model. One can see that 7 follows the same qualitative behavior as in the EMT case with 
a substantial increase upon lowering the voltage. As a conclusion, EMT and RRN models lead to very close results 
about the T c distribution and its variation with the gate voltage provided a filamentary Christmas tree structure is 
used for the RRN. 

- Multiband BCS scenario — The inhomogeneous character of LXO/STO systems, clearly apparent both in the 
normal-state [_R(T)] and in the SC [J S (T)] properties, entails a distribution of local critical temperatures with an 
average depending on the overall density (i.e., gating) T c (V g ). This dependence can be extracted from the resistance 
fits [Fig. p] and it provides insight on the the intra-puddle pairing mechanism. A recent magnetotransport analysis [1] 
demonstrated the coexistence of a rather large density of low-mobility carriers with a smaller density of high-mobility 
carriers in the SC regime of LTO/STO. In particular, superconductivity seems to be triggered by the presence of 
the high- mobility carriers. A similar behavior seems to occur in LAO/STO [21] . According to this scenario, we 
propose that the 2DEG formed at the LXO/STO interface may be described as a multi-band system [22H25] . and the 
occurrence of superconductivity may be captured by means of a BCS-like Hamiltonian (see Supplemental Material) . 
We take for the €-th sub-band of our 2DEG a parabolic dispersion law 

2me lX 2mt, v ' 

where m is the (possibly anisotropic) effective mass of the charge carriers. Following the evidence that superconduc- 
tivity is related to the appearance of high-mobility carriers, we represent the whole set of low-lying bands with one 
sub-band [£ — 1) collecting all the low-mobility carriers with vanishingly small T c , while the high- mobility carriers in 
the sub-band £ = 2 give rise to a finite T c . Thus, SC puddles are regions where the sub-band £ — 2 is locally filled, 
whereas the (weakly localizing) metallic background corresponds to regions where the sub-band £ = 2 is empty. To 
reproduce the experimental results, we are led to a suitable choice of the pairing amplitudes gu> (intraband for £ — £' , 
interband, for £ ^ £'): gn <C (312,321) <C 322 (this condition is consistent with the analysis of a two-band model in 
Ref. l26|) . According to the standard BCS approach, the pairing amplitudes are only effective in a limited dynamical 
range |ek,£ — A*| j |£kV — /-* | < hcoo, where w is a characteristic frequency of the mode that mediates pairing and /1 is the 
chemical potential. We assume henceforth that the bottoms of the two sub-bands are well separated, £2 — £1 ^> fi^o, 
and take £2 = as the reference energy level. 

While more general (though standard) expressions are derived in the Supplemental Material, here we only report 
the expressions obtained for 311 = 312 = 321 = 0, so that the system is not SC until the carrier density reaches the 
value such that fi = 0. For < /1 < haj , T c increases with increasing fi 

T c ^l.U^/huj Q fie- 1/X2 , 

where we merged into the single dimensionless coupling A2 = N^§22 the intraband 322 coupling and the density of 
states N$ of the £ = 2 sub-band. For fi > hu j , T c s aturates to its maximum (BCS) value T™ ax w 1.1A huj e' 1 ^ 2 . 
Thus, T c (ji) = 0, for n < 0, T c (/z) = T™ ax yj[i/hw , for < [i < hu , T c (p) = T™ ax ', for n > hw , and the range 
of variation of /1 which corresponds to an increasing T c is a direct measure of the characteristic energy scale of the 
pairing mediator, huiQ. 

We now proceed to fit the curve T c (V g ) extracted from the experimental data by means of EMT, with our theoretical 
curve T c (fi). The relation between V g and fj, is only approximately linear and its determination is described in detail 
in the Supplemental Material. The resulting critical temperature for various values of V g is reported in Fig. p5] (blue 
dashed line and crosses). From the fit, we finally obtain the dimensionless coupling constant A2 ~ 0.25 (which, 
consistently with our assumption, is in the weak coupling regime) and Hloq ~ 23 meV, which is the energy typical of 
phonons in STO[27], 
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FIG. 3. T c extracted from the EMT fit to experimental data (red solid line and empty circles) and the theoretical T c calculated 
within the two band model (see text), as a function of the gate voltage V g (blue dashed line and crosses). 

In conclusion, wc describe the superconductivity in LAO/STO and LTO/STO oxide interfaces within a scenario 
in which SC puddles form a percolating network. In this framework we showed here that the sheet resistance in 
LTO/STO interfaces is very well described by EMT or by a Christmas tree model of RRN for an inhomogeneous 
2DEG with a substantial filamentary character. Fitting the experiments allows to extract the random distribution 
of T c at various carrier densities (or V g ). A similar approach is also adopted to fit the micrometrically averaged 
superfluid density in LAO/STO |S]. Assuming an effective two-band model with superconductivity triggered by the 
presence of few high-mobility carriers in a higher band, we account for the density dependence of the intra-puddle T c 
within a simple BCS weak coupling scheme. As an important by-product we find that the range of variation in V g of 
the (average) intra-puddle T c is directly related (via /i) to the typical energy of the pairing mediator ojq. 

S.C. and M.G. acknowledge financial support from "University Research Project" of the "Sapienza" University 
n. C26A1I5HTN. This work has been supported by the Region Ile-de-France in the framework of CNano IdF and 
Sesame programs, and by the DGA PhD program. Research in India was funded by the Department of Information 
Technology, Government of India. 
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Supplemental Material for Inhomogeneity and multiband 
superconductivity at the LaTiOs/SrTiOs and LaAlOs/SrTiOs interfaces 

In this supplemental material we provide further detail about the following two aspects of our work: the random 
resistor network used to model the inhomogeneous character of the oxide interfaces and to fit the experimental 
resistivity curves (section A) and the two-band model used to describe superconductivity inside the superconducting 
(SC) puddles (section B). 

A. Details on the Christmas-tree model 

To describe the inhomogeneous electric conductivity of LaTiOa/SrTiOa and LaAlOa/SrTiOs interfaces, we solve a 
random resistor network (RRN) consisting of two types of bonds: (i) normal bonds with a constant finite resistance 
Ro = 1 and (ii) SC bonds with a temperature-dependent resistance R l = Rq9(T — T C J ), where T* is the critical 
temperature of bond i and 0(x) is the Heaviside step function. The critical temperatures are random variables obeying 
a Gaussian distribution. The SC bonds form a spatially correlated cluster generated through a growth process known 
as diffusion limited aggregation (DLA) on which patches of circular shape with a diameter about 30 bonds long are 
superimposed, leading to a structure as shown in figure lf\ (for details on the cluster growth see Ref. ITS]) . The number 
of patches is chosen to yield RNNs with fractions of SC to total bonds equal to w — 0.3, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7. 

Solving numerically Ohm's law on the bonds and Kirchoff 's law for current conservation at the nodes of the network, 
the resistivity curves R W (T) are obtained. For a better statistic, we calculate R W (T) for 5 different realizations for each 
weight w and take the arithmetic average. These curves are then used to fit the experimental curves Ry measured 
at voltages V = 50, 70, 100, 150, 200 [V] in the following way: The high-temperature sheet resistance of the numerical 
solution is rescaled to the experimental one, with the exception of 100[V], 150[V] and 200[V], where, to avoid the high 
temperature step, the fit is only done for the temperature interval [0.03,0.21] and the sheet resistance is rescaled to 
Rv(T — 0.21). The numerical temperature is then redefined to T = (T — a)/b, where a and b are chosen so as to 
minimize the difference 

d = Y J \Rv{T,)-R w {T j )l (3) 

the sum running over all measured temperatures. The rescaled curve which minimizes the difference d, R w =w{T) 1 
is taken as the "correct" one and the sample at the corresponding voltage is assumed to have a weight of SC bonds 
equal to w. 

In addition to the quasi ID long-range connectivity and quasi 2D short-range connectivity, the Christmas-tree model 
has the characteristic of closed loops covering almost entirely the cluster. With the size of a bond about 10 nm, these 
loops are of the order of a micrometer. This aspect is relevant when interpreting superfluid density measurements. 
Indeed, the question arises about the magnetic response measured when performing a SQUID experiment, the SQUID 
pick-up coil recording the local magnetic field produced by a field coil that extends over a few micrometers. Under these 
conditions, the measurement resolution exceeds the (average) size of a loop and the measurement cannot distinguish 
whether a closed loop encompasses a SC or a non-SC region. More precisely, the SC loops being smaller than the 
experimental resolution, the SQUID sees all bonds inside a loop as if they were SC. As a result, the SQUID experiment 
can yield a superfluid fraction w ~ 1 all the while the transport measurements yield resistivity curves that can be 
fitted with w ~ 0.5. 

While fitting the experimental data may seem to require simply a cluster with a small long-range connectivity and 
a high density, we stress that the issue is more subtle than that. Indeed, we considered a series of clusters sharing 
some characteristics of the Christmas-tree model and calculated the corresponding resistivity curves (see figure |8| . 
Starting from the assumption that the cluster should be filamentary, dense and possibly exhibit loops, we considered 
a symmetric random walk (SRW) as shown in figure El Even though the connectivity is rather small, it is too large 
to yield pronounced tails as experimentally observed. In order to make the system more ID and possibly denser, we 
reduced the number of random walks and superimposed SC patches (obtaining a cluster similar to the Christmas-tree) 
as shown in figure [3J One does obtain a tailish feature at low temperature, which however is not in agreement with 
the experimental data. Moving on to a DLA cluster without patches (figure [6]) , we found resistivity curves leading 
to rather satisfactory fits of the transport data. However, this cluster does not exhibit loops and, as a theoretical 
drawback, one cannot easily tune the density with keeping the long-range connectivity low. Finally, when adding 





FIG. 4. Symmetrized random walk cluster without patches;FIG. 5. Symmetrized random walk cluster with patches; weight 
weight w — 0.2. 



w = 0.5. 
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FIG. 6. Diffusion limited aggregation cluster without patches;FIG. 7. Diffusion limited aggregation cluster with patches; 
weight w ~ 0.3. weight w = 0.6. 

patches, these two limitations are overcome, which leads us to conclude that our Christmas-tree model captures the 
main features of the SC state at the interface. 



B. Two-band superconductivity for LXO/STO interfaces 

We propose that the 2DEG formed at the LTO/STO interface may be described as a multiband system, and the 
occurrence of superconductivity may be captured by means of a BCS-like Hamiltonian 






1 V^ t 



k,£,t C -M,i C -k',£',i C k',£',t 



(4) 
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Temperature 

FIG. 8. Fit of resistivity curve measured at V — 100 [V]. Grey curves report the experimental data Rv=ioo (duplicated and 
shifted for ease of reading), coloured curves report (from left to right) numerical calculations obtained for DLA with patches, 
DLA without patches, SRW with patches, SRW without patches. 



where c^ t CT (c k £ a ) creates (annihilates) an electron with two-dimensional wave vector k, parallel to the plane of the 
interface, and spin projection a, belonging to the ^-th sub-band, with dispersion law 



£k,. 



£e + 






2m 



e,x 



2mi, 



Hereafter, eg label the levels arising both from the level quantization inside the confining potential well and from the 
various bands mostly due to the ti a orbitals of Ti in STO. m is the (possibly anisotropic) effective mass of the 2D 
charge carriers, fi is the chemical potential, N is the number of k vectors allowed by the boundary conditions in the 
first Brillouin zone of the two-dimensional lattice. The parameters gu' are the (intraband for £ — £', interband for 
£ 7^ £') pairing amplitudes, and the tilde in Eq. Q indicates that the sums are restricted to electron states such that 
\eu.,£ ~ A 4 1) kk',£' — lA < hujo, where loq is a characteristic frequency of the mode that mediates pairing. 

In our description of the LTO/STO interface, we assume, for simplicity, that only two sub-bands, labeled by £ = 1, 2, 
are involved. The £ = 1 sub-band contains the low-mobility carriers found in Refs. H]and[T2 whereas the £ = 2 sub- 
band accommodates the high-mobility carriers. The SC puddles are regions where the sub-band £ = 2 is locally filled, 
whereas the metallic (or weakly localized) background corresponds to regions where the sub-band £ = 2 is empty. To 
reproduce the experimental results, we are led to assume that gu <§C (712,521 *C 322 (this condition is consistent with 
the analysis of two-band model in Ref. [551 but different assumptions gu < (712,321 or 312,321 ^ 322, would lead to 
similar conclusions as far as the present analysis is concerned). Finite 312,321 guarantee that superconductivity is 
induced in the first sub-band, as soon as it establishes in the second sub-band. 

The equation that determines the SC critical temperature T c in a two-band system is 



[1 - 3iiIIi(T c )] [1 - 32 2 n 2 (T c )] = 3i232ini(T c )n 2 (T c ) 



where 



n*(T) 



J ^(z)tanh(|) 



ji = (kbT)^ 1 , Ni(z) = N^/d(z — fisi + fin) is the density of states of the £-th sub-band (in the dimensionless variable 
z = (3e), with N° = a 2 ^/m,£^ x m£^ y /(2Trh 2 ), a being the spacing of the two-dimensional lattice. There are four possible 
regimes for IT(T): if ei — fj, > Huio, the chemical potential falls below the bottom of the band and there are no states 
available for pairing within the shell of width 2hu>o, hence IT(T) = 0;if0<££ — /i< hujQ, the chemical potential still 
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falls below the bottom of the band, but there are some states available for pairing within the shell of width 2Hu>q, 
hence IT(T) = N° hi[y/hu)o/(ei — //)]; if — Hu>o < ££ — (J, < 0, the chemical potential falls within the band, although 
there are some states unavailable for paring within the shell of width 2hu} , hence II^(T) = N® ln[1.14/3y (/i — et)tiUo]; 
if &i — fj, < —hujQ, the chemical potential falls within the band and all the states within the shell of width 2hojQ are 
available for pairing. In this case the standard BCS result He(T) = N® \n(l. 14 f3huJo) is recovered. We do not report 
here the crossover expressions when \et — /x| < kbT; we also assumed that fiwo 3> kbT, since we shall show that our 
system is indeed in the weak-coupling regime (i.e., N^gw "C 1 for all the coupling constants gw). 

We assume henceforth that the bottoms of the two sub-bands are well separated, e 2 ~ ?i > ^Woj an d take £2=0 
as the reference energy level. For gu = 0, the system is not SC until the filling reaches the value such that /i = — Tuvq. 
However, T c will stay exponentially small, until jjl = 0. In particular, letting \x = —Huq + Sfi, for Sfi — > + we find 



T c « 1.14 fkj e 



-2hui /(N°N°g 12 gii&[i) 



In the following, we assume that T c — in this regime. In the range < /i < Hluq, T c becomes sizable and increases 
with increasing /z. Neglecting small corrections due to 512 and 321 one finds 

T c « 1.14 y/KJ^ic- 1 /^ 9 ^, (5) 

while, for \i > hujQ, T c saturates to its maximum (BCS) value 

T™ ax «1.14&Joe _1/(JV2 ° ff22) . (6) 



Thus, T c {n) = 0, for fi < 0, T c (p) = T c nax ^J h/Huq, for < \x < hio , T c (p) = T™ ax , for /1 > fiw , and the range 
of variation of /1 which corresponds to an increasing T c is a direct measure of the characteristic energy scale of the 
pairing mediator, Itioq. 

The chemical potential /1 is related to the variation of the carrier density 8n in the second sub-band as n = Sn/ftN®). 
On the other hand, we extract the value of Sn(V g ) from a self-consistent calculation of the sub-band structure in the 
quantum confining potential [4l, in the presence of the gate voltage V g . The procedure is the following: for each V g 
the full numerical solution of the coupled Schrodinger and Poisson equations yields the sub-band structure in the 
confining potential; the Fermi level Ep is then fixed to accommodate all the carriers; the SC sub-band is identified as 
the one that starts to be filled in correspondence of the value of V g at which superconductivity is observed to occur; 
the chemical potential /i entering in our two band model is obtained as the difference between Ep and the bottom of 
the SC sub-band. 

We are thus able to extract the dependence /J,(V g ) to describe the chemical potential entering the £ = 2 band of the 
SC carriers and obtain the fit to the experimental data shown in Fig. 3, for the set of parameters considered in Ref. 
[3], where the sub-bands have all d xy character, with light isotropic mass m XtV « 0.7 mo . 



